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Non-Markovian quantum state difTusion (NMQSD) provides a powerful approach to the dynamics 
of an open quantum system in bosonic environments. Here we develop an NMQSD method to study 
the open quantum system in fermionic environments. This was a long-standing problem involving 
anticommutative noise functions (i.e., Grassmann variables) that are intrinsically different from the 
noise functions of bosonic baths. We obtain the NMQSD equation for quantum states of the system 
and the non-Markovian master equation. Moreover, we apply this NMQSD method to single and 
double quantum-dot systems. 



Introduction. — The theory of open quantum systems 
has become an increasingly important topic in, e.g., 
quantum information science, quantum measurement, 
and quantum optics. Traditionally, the dynamics of an 
open quantum system was often investigated using a 
Markov master equation derived by invoking the Born- 
Markov appoximation. However, this formalism fails for 
many solid-state systems (see, e.g., [l|) where the system- 
environment coupling is strong and the environment is 
structured. Thus a non-Markovian master equation is 
required when considering the memory effect and back 
action of the environment. It is known that the derivation 
of an exact non-Markovian master equation has long been 
a challenging task. One of the breakthroughs is the ex- 
act non-Markovian master equation for quantum Brown 
motion model derived by Hu et al. [4I using the Feynman- 
Vcrnon influence functional path-integral method [3|. 

Of all the theoretical strategies used to deal with open 
quantum systems, a non-Markovian quantum trajectory 
theory known as non-Markovian quantum state diffusion 
(NMQSD) [J,|5| provides a powerful approach to the dy- 
namics of an open quantum system in bosonic environ- 
ments. In this approach, when the so-called O-opcrator 
is obtained, the quantum dynamics of an open system 
is determined by solving the NMQSD equation (i.e., a 
diffusive stochastic Schrodingcr equation) and the non- 
Markovian master equation can also be derived [6|, lZ[ . In 
contrast to the conventional master equation under the 
Born approximation, this non-Markovian mater equation 
is derived non-pcrturbatively, so it applies even for a 
strong system-environment coupling. Indeed, some ex- 
act O-operators have been found in a variety of quantum 
models [j49i, including multilevel models [10|. 

In addition to bosonic baths, fermionic baths are also 
involved in many physical systems, particularly in solid- 
state systems. The Feynman- Vernon influence functional 
path-integral method can also be used to study the quan- 
tum dynamics of an open system in a fermionic environ- 
ment (see, e.g., [11|). Recently, this path- integral method 



tion of a doublc-quantum-dot system [l^. Nevertheless, 
the extension of NMQSD method to an open quantum 
system in fermionic baths has been a long-standing un- 
solved problem because this open system involves anti- 
commutative noise functions (i.e., Grassmann variables) 
that are intrinsically different from the noise functions 
of bosonic baths. In this Letter, we develop an NMQSD 
method to study the open quantum system in fermionic 
baths. This NMQSD approach is formulated in a non- 
perturbative manner and it applies for both weak and 
strong system-environment couplings. We not only ob- 
tain the NMQSD equation for quantum states of the 
system, but also derive the non-Markovian master equa- 
tion. Moreover, as interesting examples, we apply this 
NMQSD method to single and double quantum-dot sys- 
tems. 

Quantum state diffusion equation. — We consider a 
quantum system coupled to two fermionic baths: H = 



H, 



H, 



Hint, with (wc set h ~ 1) 



was extended to derive the non-Markovian master equa- 



k 

Hint = '^{gLkclaLk + gRkC^ajik + H.c). (2) 

fc 

Here iJgys denotes the Hamiltonian of the system, -ffenv 
is the Hamiltonian of the two electric leads acting as 
fermionic baths, and Hint niodels the interactions be- 
tween the system and the two baths. The spectral den- 
sity function of each bath is J\{uj) = ^i.\gxk\^S{uj—ijj\k), 
where X = L or R. In Eq. ([1]), a^j, {a\k) is the fermionic 
creation (annihilation) operator for a quantum state with 
wave vector k in the left or right lead. We assume that 
the system of interest couples to the two leads via single 
channels characterized by the fermionic creation (annihi- 
lation) operators c^j, (cxk) [see Eq. ^]. Extension to a 
multi-channel case is straightforward. 

In an NMQSD approach, environments are required to 
be initially at zero temperature, so as to conveniently 
represent the environmental degrees of freedom with 
the coherent state basis. As for environments initially 



with a nonzero temperature, one can map the nonzero- 
temperature density operator to a zero-temperature den- 
sity operator using a Bogohubov transformation |13| . 
In the case of fermionic baths, this requires to add 
^\k'^>^kb\kb\j; to Eq. ([1]), corresponding to the part in- 
volving holes in the electric leads. The Bogoliubov trans- 
formation for fermionic operators can be introduced as 



axk = Vl - nxkdxk 
b\k = Vl - nxkexk - 



/nxke'^^ 



(3) 



where nxk = [e'^'^^'=~''^)/'^^^ + l]^^ is the average number 
of electrons in the fcth state of the left (right) electric lead 
with chemical potential fix- In Eq- &, the coefficients 
are determined by the requirement that the derived mas- 
ter equation reduces to a Lindblad form in the Markovian 
limit. The transformed Hamiltonian Ji is written as 

n = H^ys + Y^ \ojxk{dlkdxk + exkeik) 

Xk 

+ {V^9xk(^\exk + Vl - nxkgxkc\dxk + H.c.) , (4) 

where the new fermionic operator dxk i^xk) corresponds 
to the annihilation of electrons (holes) in the virtual 
fermionic baths. Note that the effects of temperature 
are incorporated into the transformed Hamiltonian and 
the fermionic baths with nonzero initial temperatures are 
mapped to virtual fermionic baths with zero initial tem- 
perature. 

In the interaction picture with respect to the environ- 
mental Hamiltonian "Hcnv = Hxk^>^k[d\j,dxk + eAfce^^,), 
the total Hamiltonian reads 



1, the state |^t) can be expressed as 



l*t> = e 



-z z — w w 



zw) <»\'ipt{z*,w*))d^zd^w, (8) 



where z*z = Y^xk^xk^xk, w*w = ExkKk^xk, d^z = 
Y[xk^-^xk'^^>^k, and d'^w = Y[xk'^''^xk^''^>^k- The actions 
of annihilation (creation) operators dxk and exk {dx^. 
and el, J on fermionic coherent states satisfy the rela- 



-^Xk/ 



tions [T3|; dxk\z)x = zxk\z)x, dljz) 
exk\w)x = ■wxk\w)x, and e|^|w)A 



AfcFM 



dz\k 



\z)x, 

^Afci-/A - -&;;j7rk>A- When 
projecting onto the coherent-state basis, the equation of 
motion (|6]) can be reduced to the NMQSD equation for 
a pure state of the system \tpt{z*,w*)) = (zwl'it)- 



|l^*> 



iH^ysH't) - Y. [cxzlitMt) + c{wl{t)\ijt 



V4 / aAi(i-s) 
^ Jo 



Szlis) 



V'ca / ax2{t~ s) 

^ Jo 



Swlis) 



\ipt)ds 



\ipt)ds, 



(9) 



which initiates from \ipt=Q{z* ,w*)) = \fo)- Here the 
noise functions z^(t) and w^{t) are defined as 



k 



z*xit) = -*y" V^l - nx{uJxk)9*xkZ*xke" 



lit) = -i^ v^ 



uJxk)9xkWxke 



(10) 



T-L{t) = Hsys + V(V^^Afc.gIfcCAeAfce 



icJxkt 



Xk 



+V1 - nxkgxkAdxke^'^^"' + H.c), (5) 

and the quantum state of the total system satisfies the 
equation of motion 



dt\^t)=-i'H{t)\^t 



(6) 



We assume that the quantum state of the total system 
is factorized at the initial time i = 0, so that j^^o) ~ 
Iv'o) <8) |0), with the virtual fermionic baths initially in the 
ground state (i.e., at zero temperature): |0) = {^;)j|0)Ad'8' 
|0)Ae, where dxk\G) = 0, and eAfc|0) = 0. 

Define a femionic coherent-state basis for the environ- 
mental degrees of freedom: 



\z)x 



\zk)x 



\Wk)x 



zw) = {8);^|z)a ® \w)x, with 



-^Zftii-Afceif 



|0), 



(7) 



where Zk and Wk are Grassmann variables that obey the 
anticommutation relation. With the completeness rela- 
tion for coherent states / e~^ ^~™ '^\z'w){z'w\d'^zd'^w = 



The temperature-dependent environment correlation 
functions are axi[t — s) = M{zx{t)zl^{s)} = 
j'du[l - nA(w)]JA(w)e-^"(*-^), and ax2(t - s) = 
M{wx(t)wl{s)} =/dwnA(a;)JA(w)e*'^(*-"), where X{-} 
denotes the statistical mean over all noise variables: 
M{-) = /e-^*^-™'™{-}d2zd2u,. 



Introducing 0-operators by 



:\Mz\w*)) 



Oxi{t,s,z*,w*)\Mz*,w*)), and j:;^\^t{z\w*)) = 

Ox2{t,s,z\w*)\il}t{z\w*)), we can write the NMQSD 
equation in a time-local form: 



d_ 
di 



IV't) - ~iH,ys\i^t) - Y h4(i) + cWxit) 
X 

+cldxiit,z*,w*) + cxOx2{t,z*,w*)\ \^Pt), (11) 



where Oxn = /g dsaxnit - s)Oxn{t,s, z*,w*), n == 1,2. 



With the consistency conditions 



d JWh 



d\-i>t 



d &\'4't) _ S d\'ipt) 
dt Sz'^{s) 5z^{s) dt 

as well as the initial condi- 



tions Oxi{t, s, z* ,w*)\t=s= cx and Ox2{t,s,z*,w*)\t=s= 
c\, we obtain the equations of motion for the O- 



operators: 
dt 



A' 
A' 

(12) 



where the square and curly brackets denote the com- 
mutator and anticommutator, respectively, and Qn = 






^t SOri 
-R (5A„ 



^ <5A„ 



___ c^fxf ' ^ith Ai = z*(s) 

and A2 — wj(s). 

Master equation. — The reduced density operator of an 
open quantum system by tracing over the environmental 
degrees of freedom can be obtained by taking the sta- 
tistical mean on a density operator related to the state 



pt = Tre„v|*t>(^t|=7W{Pt}, 



(13) 



where Pt = lijjtiz* ,uj*)){tpt{—z,—w)\. Using the relation 

and Eq. pT|) . we derive the following non-Markovian 
master equation: 



^ = -z[H,ys, Pt] + E {[cx,M{PtOUt, -z, ~w)}] 



-[cx,M{Ox2it,z*,w*)Pt}] 
+ [clM{PtOUt,~z,-w)}] 



(14) 



This master equation is derived non-perturbatively, so it 
applies even for a strong coupling between the system 
and the environments. Moreover, in addition to trace 
preserving, it also preserves the positivity and hcrmitic- 
ity. 

In the Markovian limit, axi{t — s) — > [l^hx\T\6{t~ s), 
and ax2{t - s) -^ fixTxS{t - s), where Tx = 2'!rpx\gx\'^, 
with X = L [R], is the electron tunneling rate between 
the system and the left (right) lead. Also, the time- 
integrated 0-operators become Oxi -^ ^^^(l — nA)cA, 

~ 1 — + 

and 0x2 -^ ^^x^xCx- Therefore, the master equation 

Til) is reduced to 



dpt 
dt 



Tx 



+ E -^ nx{2c\ptcx - cxc[pt - Ptcxc[) 

A 

+ (1 - nx){2cxptcl - c\cxpt - Ptc\cx) ■ (15) 

It is clear that this Markov master equation has a Lind- 
blad form. 

Application to quantum-dot systems. — Below we apply 
our NMQSD approach to single and double quantun-dot 
systems. 



(i) Single quantum dot. Suppose that the single quan- 
tum dot is in the strong Coulomb blockade regime, so 
that only one electron is allowed therein. The Hamil- 
tonian of the system is written as -ffsys = ^qc'^c, and 
cl — cn = c for iJjnt in Eq. ^. The non-Markovian 
master equation is exactly derived as 



^=-z[i7sys,Pt]+ri(t)[c,ptct]+r2(t)[c,cVt] 



rut)[c\cpt]-r*2{t)[c^,ptcl 



with time-dependent rates 



(16) 



Tj{t)^ / [ai{s-t)Aj{t,s)-a2it-s)Bj{t,s)}ds, (17) 
"'0 

where aj{t) = aLj{t) + afi,j{t); Aj(t,s) and Bj{t,s) are 
determined by the integro-differential equations: 

{—-ioJo)Aj{t,s)+ f P{s-s')Aj{t,s')ds' = U{t,s) 
OS Jo 

(— - iuJo)Bj{t,s) + J /3(s - s')Bjit,s')ds' = Vit,s), 

(18) 

with /3(s— s') = ai(s' — s)-f 02(5 — s'): U{t, s) = J^ [02(5 — 
s')h{t,s')ds\ V{t,s) = /o[ai(s' - s)h{t,s')ds\ and the 
final conditions at s = i: Ai{t,t) = Bi{t,t) = 1, and 
^2(i, t) = B2{t, t) = 0. Here h{t, s) satisfies the equation 
(^ - iujo)h{t,s) ~ X*^(s - s')h{t,s')ds' = 0, with the 
final condition h{t,t) ~ 1. 

In the Markovian limit, ri(i) ^ i[l - nL{u!o)]TL + 
|[1 - nR{^o)]TR, and r2(i) -^ -~^nL{uJo)TL - 
^fifi{uJo)Tf(. Let us consider the zero-temperature case 
with 71^(1^0) ^^ ^(ma ~ '^o); were is the Heaviside step 
function. If the single-dot level uiq lies within the energy 
window fiL > luq > fifj, the master equation (|16p reduces 
*° SiPt "^ -iuJo[c'^c,pt] + \TL{2c^PtC- cc^pt - Ptcc^) + 
\Tii{2cptC^-c^cpt-ptC^c). With the basis state 1 0) (|1)) 
which denotes an empty (occupied) dot, it follows that 
the density matrix elements satisfy 



Poo = — ElPoo + r^pii, 

Pii ~ r^poo — Tjipii, 

Pio = -(iwo + Tl + rfl)pio. 



(19) 



which are exactly the rate equations obtained by Gurvitz 
and Prager [15j . 

(ii) Double quantum dot (DQD). Suppose that the 
DQD is in the strong Coulomb regime, so that at most 
one electron is allowed in each dot. The Hamiltonian of 
the DQD can be written as i^sys = ujic[ci + UJ2C2C2 + 
rio(c2Ci + c[c2), where fto denotes the inderdot coupling. 
For Hint in Eq. ^, cl ~ ci, and cr = €2- The exact 
non-Markovian master equation is given by 



dpt 
dt 



■ - i[Hsys,pt] +TLi{t)[ci,ptc\] +TL2{t)[ci,c\pt] +TLz{t)[ci,ptC2] +TLi{t)[ci,c\pt] 

+ TRl{t)[c2,ptc\]+TR2{t)[c2,c\pt]+TR,3{t)[c2,ptcl]+TBA{t)[c2,c\pt\ 

-r2i(t)[c|,cipt] ~Tl2(t)[c\,PtCi] -Tl^{t)[cX,C2Pt] -Tl^{t)[c\,ptC2] 

- ^*Rl{t)[c\,CiPt\ - T*R2{t)[cl, PtCi] - r^3(i)[4>C2Pt] - r|j4(t)[4,ptC2], 



(20) 



r 



with tiiTLC-dcpcndcnt coefficients 



rAj(i) 



[a\i{s-t)Axj[t,s) - ax2{t- s)Bxj{t,s)]ds, 

(21) 



where A\j(t,s) and B\At,s) satisfy a set of integro- 
differential equations [l6|, with the final conditions: 
Aii(i,i) = Afl3(i,i) = BL2{t,t) = BRi{t,t) = 1, and 
Axj{t,t) = 0, Bxj{t,t) = for other A and j. A similar 
non-Markovian master equation was also obtained using 
the Feynman- Vernon infiuence functional path-integral 
method [H. 

In the Markovian limit, TLi{t) — >■ ^[1 — nL{uJi)]TL, 
Tri -> 0; ri2(i) ^ -iflL(wi)rL, rR2 -^ O; TLsit) -^ 
0, Tfls ^ ^[1 - nR{LU2)]rR; and TLi{t) -> 0, Tra -^ 
— ifi,fl(a;2)]rfl. We also consider the zero-temperature 
case with hx{i-^n) ~^ ^(/^A ~ ^n)-, a-nd the two single-dot 
levels of the DQD all lie within the energy window p^ > 
Wn > PR, where n = 1,2. We use |/), Z = 0,1,2 and 3, 
to denote the states with both dots empty, the left dot 
occupied, the right dot occupied, and both dots occupied, 
respectively. From Eq. ([20)1 . it follows that the master 
equations for density matrix elements are reduced to 

Poo ~ — TlPOO + r_RP22, P33 = rLP22 — ^RP33, 

Pll = TlPoO + r_RP33 + iOo(/3i2 - P2l), 

P22 = -(Fl + Tr)p22 - i^o{Pl2 - P2l), (22) 

Tl + Tr 



Pl2 



-i{uJi - uj2)pi2 + «^o(pii - P22) 



-Pl2, 



which are identical to the rate equations obtained in |12l . 
15[. For a DQD, both intradot and interdot Coulomb 



repulsions can play an im por tant role in the Coulomb- 
blockade effect (see, e.g., [17|). Thus, if both intradot 
and interdot Coulomb repulsions are so strong that only 
one electron is allowed in the whole DQD, the master 
equations for density matrix elements are reduced to 

Poo = — r^poo + ^RP22 

Pll = TlPoo + i^o{pi2 - P21) 

P22 = -rflP22 - i^o{pi2 - P21) (23) 

P12 = -«(^i - ^2)Pi2 + i^aiPii - P22) ^/'i2' 

which are exactly the rate equations obtained in |18l |. 



Conclusion. — Wc have developed an NMQSD method 
to study the dynamics of an open quantum system in 
fermionic baths. We not only obtain the NMQSD equa- 
tion for quantum states of the system, but also derive the 
non-Markovian master equation. This non-Markovian 
approach is formulated in a non-perturbative manner and 
it applies even for a strong coupling between the system 
and the fermionic baths. Moreover, as useful examples, 
we have applied this NMQSD method to single and dou- 
ble quantum-dot systems. 
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